Photoisomerization in a Dissipative Environment 
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We investigate photoisomerization (PI), the shape change of a molecule upon photoabsorption, 
in a dissipative environment using a simple spin-boson model. We identify two classes of environ- 
ment depending on whether it "entangles" with the molecule. In the absence of entanglement, the 
environment merely causes a blue shift of the required photon frequency and reduces the quantum 
efficiency of PI. With entanglement the molecule can undergo a quantum phase transition between 
a state that photoisomerizes to a state that does not. 

PACS numbers: 85.65.+h, 33.80.-b, 42.88.+h, 82.30.Qt 



Photoisomerization (PI) is a phenomenon where a 
molecule absorbs a photon and changes its shape 
(conformation). 0, 0, H EL H 

This important phe- 
nomenon is relevant to vision and other biological pro- 
cesses. Recently it has been proposed that molecules 
which photoisomerize (the "PI molecules") can be used 
as light sensitive molecular switches [1, 0| and optome- 
chanical devices. Q A good example of PI molecule is 
"azobenzene" . [l| It switches between two isomers, the 
trans and cis configurations upon absorbing light. 

For microscopic systems, environment usually plays a 
critical role in affecting their quantum mechanical behav- 
ior. PI of molecules is not an exception. The goal of this 
paper is to study PI in the presence of an environment. 
Although the rest of the discussions are coined in terms 
of PI, our theory should be applicable to a variety of sit- 
uations. Examples include Jahn- Teller distortion in an 
environment and quantum bit in an environment. 

The physical mechanism of PI has attracted a lot of in- 
terests in the past two decades. In this work we present a 
simple theory for it. Our theory differs from those in the 
literature in that we emphasize the importance of the en- 
vironment. Specifically, we suggest the existence of two 
types of environment, depending on whether it "entan- 
gles" with the PI molecule. Without entanglement the 
environment merely causes a blue shift in the photon fre- 
quency necessary to induce PI, and reduces the quantum 
efficiency. With entanglement the molecule can undergo 
a quantum phase transition between a state that photoi- 
somerizes to a state that does not. Our investigation of 
the environmental effects on PI is stimulated by a very 
recent scan tunneling microscopy experiment which re- 
ported the loss of PI functionality when azobenzene is 
in close contact with a metallic (gold) surface, and the 
regaining of the PI functionality as the molecules are el- 
evated sufficiently far from the surface. 

The irreducible in gred ients of a PI molecule are (1) 
two electronic states[10( and (2) a structure degree of 
freedom, and (3) a strong coupling between the two. We 



model (1) by a quantum bit (or a quantum spin) where 
the up and down states correspond to the two electronic 
states in question, and (2) by a simple harmonic oscil- 
lator. The coupling between the two sets of freedom is 
described by the following Hamiltonian 
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Here x and p are the dynamic variables of the harmonic 
oscillator, a z is the third Pauli matrix, and 2to describes 
the energy difference between the two electronic states. It 
is instructive to first analyze Eq. ([1]) in the limit m — > oo 
where the harmonic oscillator is classical. In this limit, 
depending on the value of a z (±1), the energy landscape 
for the structural degree of freedom (U(x)) takes on dif- 
ferent shapes. This is shown as the red and green curves 
in Fig. 1(a). In conventional picture of PI, the system 
jumps vertically from the minimum of the red (green) 
curve to the green (red) curve upon absorbing a photon. 
Subsequently the structure relaxes until x reaches the 
minimum of the green (red) curve. The two minima of 
the green and red curves, situated at x = ±g'/k, corre- 
spond to the two relevant molecular structures. We note 
that by setting fo to zero Eq. (TTJ) describes the physics 
of Jahn- Teller effect where the structure distortion (i.e., 
x = — > x 7^ 0) spontaneously sets in to split the elec- 
tronic degeneracy and lower the total energy. 

The above intuitive picture for PI leaves unclear an 
important question: since electronic transition often pro- 
ceeds at a faster rate than structural relaxation, what 
holds the molecule from de-exciting back to the original 
electronic state during structural relaxation? We begin 
our discussions by first addressing this question. 

Restoring the full quantum mechanical nature of the 
harmonic oscillator and introducing creation (af) and an- 
nihilation (a) operators, we replace the 2nd to 4th terms 
of Eq. |T]) by <?oc z (a^ + a) + bj^a^a. This simple model 
can be solved exactly to yield the following eigen-states 



2 




a =-1 04 



< b > II 


I 1 1 1 I 1 

a c =2.5, n=o.2ci) 




A„=3, toJa B =3 - 




l x10 - 

\l. , : , 



20 40 60 80 100 



FIG. 1: (Color online) (a)The energy landscape for the struc- 
ture coordinate. (b) The optical absorption spectrum a(e) for 
the PI molecule without (red line) and with (blue line) cou- 
pling to environments. r\ is the Lorentz broadening parameter. 



and eigen-energies 
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Here a c = go/ujQ, s — ±1 (which labels the two elec- 
tronic states), and n are non- negative integers. The state 
|s,0) satisfies a\s,0) = and a z \s,0) — s\s,0). It is 
simple to show that regardless of the value for n, we 
have (±l,ro|x| ± l,n) = z \ 1 2a c implying the presence of 
two structural conformations with the associated quan- 
tum structural fluctuations. For to > 0, the ground 
state is | — 1,0) with (x) — 2a c . After absorbing a 
photon the molecule make a transition from | — 1,0) to 
| + 1, n). When that happens the structure changes from 
that corresponding to (x) = +2a c to that corresponding 
to (x) — —2a c . 

If we model the electric dipole operator as a x , the op- 
tical absorption spectrum is given by 
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Substituting Eq. © into Eq. © we obtain 
a(s) = Y^^Lofn(a c )S(e - nuj - 2t ), where 
fn(o-c) = (4ac) n e -4a <= jn\ is the Franck-Condon fac- 
tor between the initial | — 1,0) and final | + 1, n) states. 
The red curve in Fig. 1(b) shows a(e) for a specific 
choice of parameters. The closely spaced peaks originate 
from the excited states of the harmonic oscillator. The 
envelope of the red absorption curve is determined by 
the Franck-Condon factor /„(a c ). As shown in Fig. 1(b) 
the absorption practically vanishes for (e — 2t )/uj ^ 10. 
This is due to the suppression of the optical transition 
matrix element by f n (a c )- Likewise, the envelope peaks 
at the energy where the Frank-Condon factor reaches 
maximum. At this peak frequency the incident photon 
has the highest quantum efficiency in inducing isomer- 
ization. Due to the strong suppression of the optical 
matrix element, the molecule will cease to spontaneous 
emit once its energy has descended below certain energy 



(~ 10wo in Fig. 1(b)). This result suggests the presence 
of two time scales 7i,2- T\ is the time it takes for struc- 
tural relaxation to bring the molecule's energy outside 
the absorption peak, and Ti is the time for it to bring 
the molecule to the isomerized ground state. During 
t = — > T\ the molecule relaxes both structurally 
and electronically (by spontaneous emission). During 
t = T\ — > T2 the molecule can only relax structurally 
because spontaneous emission is strong suppressed. 
Thus once the photo-excited molecule survives through 
T\ it can relax to the isomerized ground state without 
the further interruption from spontaneous emission. 
This answers the earlier posed question, and at the 
same time predicts that the quantum efficiency of PI 
is determined by the probability for the molecule to 
survive through Ti.The existence of time regimes with 
different relaxation dynamics is consistent with the 
experimental results for azobenzene. 11| Moreover, due 
to the energy difference 2to, the peak photon frequency 
for the (x) — +2a c — > (x) = —2a c conversion will be 
different from the (x) = —2a c — ► (x) = +2a c one. This 
is also consistent with experiments. [ll| 

Next we consider the effects of the environment. Fol- 
lowing Calderia and Leggett[l2| we model the environ- 
ment by a continuum of harmonic oscillators. These os- 
cillators couple to the PI molecule via the Hamiltonian 
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Here a labels the environment oscillators and /1 = x or 
y or z depending on the nature of molecule-environment 
coupling. Clearly when /i = z we have [H mo \, H cnv ] = 0, 
which implies the absence of "entanglement" between the 
molecule and the environment (see later). On the other 
hand, for /i = x or y H mo i and H onv do not commute 
and the environment entangles with the molecule. As 
we shall show that these two types of environment have 
fundamentally different effects on PI. 

Let us begin with the unentangled environment (i.e. 
[L — z). In this case the eigenstates of the molecule plus 
the environment are the following direct products (hence 
no entanglement) \s, n) eg) \s, {n a }}. Here \s, n) is given in 
Eq. ^ and |s, {n a }} is given by 
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with a™ = g a /uj a . In this case the environment has no 
effect on the structure, i.e., (x), of the molecule. It merely 
modifies the absorption spectrum to 
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where D{w) = M. e - iMt+AK{ - t > , K{t) = 

/ °° ^J(o;)(e-*-l), andJH = £ a <&*(w-a;„). In the 
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case where the environmental dissipation is Ohmic[l3j we 
have J (to) = Aotoe^^^" , where Aq reflects the molecule- 
environment coupling strength, and lo c is a cut-off fre- 
quency. For, e.g, Aq = 3 and uj c /^o = 3 the absorption 
spectrum is shown as the blue curve in Fig. 1(b). It is 
blue-shifted relative to the red curve, and the peak in- 
tensity is greatly reduced reflecting a considerable reduc- 
tion of the quantum efficiency of PI. Moreover, the sharp 
phonon peaks in the red curve is now smoothened. 

Next, we consider the more interesting case where the 
environment entangles with the molecule. Such an envi- 
ronment can be simulated by putting [i = x in Eq. ([4}, 
which results in the non-commutivity between H mo \ and 
H cnv . Physically this type of environment induces tran- 
sition between the electronic states (or flips the spin). 
When resent alone, the ground state of H mo \ and H env is 



the eigenstate of a z and a x respectively. When present 
together the spin is subjected to the conflicting quan- 
tum noise exerted by the structure degree of freedom 
and the environment. If we choose the eigenstate of, 
say, a x to describe the spin, the dynamics dictated by 
-ffmoi + -Hcnv causes the spin to flip as a function of time. 
The primary effect of the environment is to induce cor- 
relation between spin-flipping events. In the following 
we use a well-established formulation [l^] which trans- 
lates above physical picture into precise mathematical 
languages. Technically we employ the imaginary-time 
path integral representation of the quantum partition 
function of the molecule plus the environment. We then 
"integrate out" the environmental degrees of freedom and 
obtain the following effective partition function for the 
molecule 
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here 2n counts the number of a x — > — a x spin flips be- 
tween r = and r = j3, the imaginary time integral 
is subjected to the contraint n < t 2 < ■ ■ ■ < r 2n , 
x(t) specifies the history of the molecular oscillator, and 
t = l/u) c . The quantity Q(tj) is an alternating sequence 
of ±1 (+1 for down — * up flip and —1 for up — > down 
flip respectively) at ti,..,t 2n - The action S[x(t), Q(rj)] 
in Eq. ([7]) is given by 



S= / dTH(r)-J2Q( 
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where V(r) = J °° duiiv~ 2 J(lo) exp(— w\t\) and H(t) = 
^[x(t) 2 + oj 2 x(t) 2 }. In the case of the Ohmic environ- 
ment (J(uj) — Aou;e^ UJ ^ UJc ) we obtain 
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In the limit where £o r o and So^o are small, a well-known 
renormalization group technique 14] has been developed 
to treat Eq. (J7J. In this treatment one progressively in- 
tegrates out pairs of up and down spin flips that are sep- 
arated farther and farther in time, and record the effects 
of such integration on the parameters of Eq. {§]). For 
V(rj — Tj) given by Eq. © the result of such calculation 
gives the following recursion relations = —2A(t 2 + 
g 2 e - % | = (1 - A/2% § = (1 - A/2)g, f = u, 
here / is the logarithm of the time separation of the 
farthest spin flip pair integrated out so far divided by 
tq, and A, t, <?, and u> are all dimensionless quantities. 



They are function of I satisfying the initial conditions 
A(0) = Ao, t(0) = t Q T , g(0) = g Q T , u(Q) = lo q t . 
These recursion relations can be solved numerically to ob- 
tain the phase diagram for the zero temperature (/3 — » oo) 
state of the PI molecule as a function of the initial io r o 
and Aq for different <?oAo an d ^>o T o- I n Fig- 2 we present 
the phase boundary between the "instanton metal" where 
the spin flips are uncorrelated (aside from the global con- 
straint of alternating up and down spin flips), and the 
"instanton insulator" where up and down spin flips are 
correlated so that they form bound pairs in time, for sev- 
eral different values of ujotq at a fixed go /to. 

To obtain information about the molecular structure 
in these two phases we compute the correlation function 
^xx{ T a — T b) = ( x ( T a) x ( T b)) ■ It can be shown that 

Txx{T a - n) = G (T a - T b ) 

+g 12 J dndriGoin - T a )G (T 2 - r 6 )C(T 1; T 2 ),(10) 

where G (t 2 -t x ) = e -^l r2 - Tl l/(2mw ) and C{t 1i t 2 ) = 
e -AS(Ti,T 2 ) w j^ n £±g the ac tion change due to the inser- 
tion of two instantons at t\ and t 2 . In the two different 
instanton (1ST) phases 
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where t s is an instanton "screening length" , and ci , c 2 , C3 
are positive constants, implying 
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FIG. 2: (Color online) The zero temperature phase diagram 
in Ao — toro plane. In constructing it we fix the value of go /to 
to 0.3. The short-dashed, dash, and solid lines are the phase 
boundaries separating the instanton metal (upper left) and 
instonton insulator (lower right) phases for uoto = 0.3, 0.6, 
and uiqTq — > oo, respectively. 



Substituting the above asymptotic behavior of C(n,T2) 
into Eq. (|10[) and noting the exponential decaying na- 
ture of Gq{ti — T2), we conclude that the molecule has 
a long/short range (x(r)x(O)) correlation function in the 
instanton metal/insulator phase, respectively. 

This suggests that in the instanton metal phase the 
molecule essentially behaves as if it is in isolation (i.e., the 
electronic eigenstates are eigenstates of a z and the struc- 
ture degrees of freedom shows non-zero (x)). In contrast, 
in the instanton insulator phase the coupling between the 
electronic state and the environment dominates. As a re- 
sult, the electronic eigenstates are essentially the eigen- 
states of a x and the two distinct molecular structures are 
averaged out. In the former case the molecule retains its 
PI functionality, while in the latter case the functional- 
ity is lost. This phase transition can be realized by, e.g., 
tuning the strength of the molecule-environment coupling 
(A in Fig. 2.). 

The transition between these two cases will be man- 
ifested as a quantum phase transition at absolute zero 
temperature. This phase transition can be detected via 
optical absorption spectrum. To compute such spectrum 
we note that in the two different instanton phases the 
spin-spin correlation function behaves 
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Since r ss (r) is the Laplace transform of the opti- 
cal absorption spectrum, i.e.,r ss (r) = J a(uj)e~' TUJ duj, 
Eq. (JTTJ) implies 



cr(w) = a 2 S(ui) + a 4 9(uj - A)(u> - A) c 



(12) 



where A = l/r s , the zero- frequency delta- function 
weight 02 is nonzero only in the instanton insulator 
phase, and 04 is a constant proportional to a\ or 03 de- 
pending on whether the system is in the instanton metal 
or insulator phase. At the quantum critical point t s 



diverges and the optical absorption gap closes. As the 
molecule loses the PI functionality a zero frequency peak 
should appear in the optical absorption spectrum. These 
features can in principle be detected optically. At non- 
zero temperatures, the quantum phase transition will be 
smeared out. In the Aq — T plane (assuming we reach 
the phase transition by tuning Aq) there will be a re- 
gion in which fc^T > l/r s (T) where the molecule exhibit 
quantum critical characteristics. 

So far we have discussed the simplest model for PI. In 
reality there could be more than one structural modes 
that couple to the electronic optical transition[15[. We 
have investigated a two-oscillator model 1(| 13] and reach 
the conclusion that so long as the two structure oscillator 
are non-degenerate, the results discussed above remain 
qualitatively unchanged. 

In conclusion we have presented a theory of photoiso- 
merization in an environment. It points to the crucial 
role the environment plays in this fascinating molecular 
quantum mechanics. 
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